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ABSTRACT 



We present a quantitative analysis of the X-ray spectrum of the echpsing high mass X-ray 
binary Vela X-1 (4U 0900 - 40) using archival data from the ASCA Solid-State Imaging 
Spectrometer. The observation covers a time interval centered on eclipse of the X-ray pulsar by 
the companion. The spectrum exhibits two distinct sets of discrete features: (1) recombination 
lines and radiative recombination continua from mostly hydrogcnic and hclium-likc species 
produced by photoionization in an extended stellar wind; and (2) fluorescent K-slioll linos 
associated with near-neutral species also present in the circumsourcc medium. These features are 
superposed on a faint continuum, which is most likely nonthermal emission from the accreting 
neutron star that is scattered into our line of sight by free electrons in the wind. 

Using a detailed spectral model that explicitly accounts for the recombination cascade 
kinetics for each of the constituent charge states, we are able to obtain a statistically acceptable 
{Xr = 0.88) fit to the observed spectrum and to derive emission measures associated with the 
individual K-shell ions of several elements. From calculations of the ionization balance using 
the photoionization code, XSTAR (Kallman & Krolik 1995), we assign ionization parameters, 
^, to several ions, and construct a differential emission measure (DEM) distribution. The 
DEM distribution spans a broad range in ^ (Alog^ ^ 2), and is centered around log^ = 2.5. 
We find that the total emission measure of the visible portion of the highly ionized wind is 
~ 3 X lO^'^ cm-3. 

The qualitative aspects of the inferred DEM distribution are consistent with a wind model 
derived from the Hatchett & McCray (1977) picture of an X-ray source immersed in a stellar 
wind with a generalized Castor, Abbott, & Klein (1975) velocity profile. Using this formalism, 
theoretical DEM distributions, parameterized only by a mass loss rate and a wind velocity 
profile, are calculated and used to predict the detailed X-ray spectrum, which is then compared 
to the ASCA data. Again, we find a statistically acceptable fit {xf- = 1-01), with a best-fit mass 
loss rate of ^ 2.7 x 10~^ Mq yr~^. This is approximately a factor of 10 lower than previous 
estimates of the mass loss rate for the Vela X-1 companion star, which have primarily been 
determined from C IV and Si IV P Cygni profiles, and X-ray absorption measurements. We 
argue that this discrepancy can be reconciled if the X-ray irradiated portion of the wind in Vela 
X-1 is structurally inhomogeneous, consisting of hundreds of cool, dense clumps embedded in a 
hotter, more ionized gas. Most of the mass is contained in the clumps, while most of the wind 
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volume (> 95%) is occupied by the highly ionized component. We show quantitatively, that this 
interpretation is also consistent with the presence of the X-ray fluorescent lines in the ASCA 
spectrum, which are produced in the cooler, clumped component. 

Subject headings: binaries: eclipsing — pulsars: individual (Vela X-1) — stars: neutron — stars: 
winds, outflows — techniques: spectroscopic — X-rays: stars 
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1. Introduction 

In wind-driven high mass X-ray binaries (HMXBs), an intense X-ray continuum fltix is produced 
through gravitational capture of a stellar wind by a neutron star or a black hole. A small fraction 
(~ 0.1%) of the total mass lost by the companion is accreted onto the compact object, and a fraction of the 
gravitational potential energy is then converted into X-radiation, which ionizes and heats the surrounding 
gas. The wind reprocesses hard X-rays from the compact object, resulting in discrete emission lines and 
continuum radiation. X-ray emission lines can be produced in a wide range of ionization states, ranging from 
cold fluorescent emission line regions to highly ionized recombination line regions. The resulting spectrum 
carries information that should, in principle, allow us to derive physical parameters that characterize the 
state and structure of the circumsource material. 

Vela X-1, the archetypal wind-fed X-ray pulsar, has been the most well-studied object in its class since 
its discovery (Chodil at al. 1967). It is an eclipsing HMXB with a pulse period of 283 s (McClintock et al. 
1976) and an orbital period of 8.96 days (Forman et al. 1973). The optical counterpart, HD 77581, is a 
BO.SIb supergiant (Brucato & Kristian 1972; Hilter, Werner, & Osmer 1972) with an inferred mass loss rate 
lying in the range (1 - 7) xlO"'' Mg yr"! (Hutchings 1976; Dupree et al. 1980; Kallman & White 1982; 
Sadakane et al. 1985; Sato et al. 1986). Its intrinsic X-ray luminosity is ^ 10^^ erg s~^, and is consistent 
with accretion from laminar flow past a stellar object (Bondi & Hoyle 1944) for a mass loss rate and velocity 
structure characterizing a typical OB supergiant. 

Prior to the launch of ASCA, X-ray emission line studios of Vela X-1 were typically limited to the iron 
K complex near 6.4 kcV; emission lines from lighter elements wore not resolved (White, Swank, & Holt 1983; 
Ohashi ct al. 1984; Lewis ct al. 1992). The solid state detectors onboard ASCA are capable of resolving 
bright spectral features from H-like and He-like ions, although, for the most part, distinguishing the effects 
produced by various emission mechanisms is still problematic. The first analysis of the ASCA spectrum of 
Vela X-1, which emphasized line emission, was presented by Nagase et al. (1994, hereafter N94). The hard 
X-ray spectrum (3-10 keV) is highly absorbed and contains bright iron K fluorescent lines. During eclipse 
(—0.1 ^ (l> ^ +0.1), the continuum radiation from the compact source is occulted by the companion, and 
emission lines that form in extended regions gain contrast with respect to the ionizing continuum, resulting 
in the line-dominated appearance of the soft X-ray band (0.5-3 keV). N94 tentatively identified most of 
the bright emission lines, while some of the fainter mid-Z fluorescent lines were identified in a later article 
(Nagase 1997). 
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Although it was suggested in N94 that the more highly ionized X-ray emission lines observed in the 
Vela X-1 spectrum are formed by cascades following radiative recombination, those authors did not have 
access to an appropriate, self-consistent, recombination spectral model suitable for a detailed quantitative 
analysis of the emission line intensities. While it is all but certain that many X-ray lines in HMXBs are 
recombination lines, it is, in fact, difficult to distinguish recombination from other emission mechanisms 
in CCD spectra, thereby disallowing a quantitative interpretation of the spectra. Liedahl & Paerels 
(1996, hereafter LP96) applied a suitable rcconibination cascade model to the interpretation of the ASCA 
spectrum of the HMXB, Cyg X-3. Infrared line spectroscopy of this object indicates the presence of a 
dense stellar wind from a Wolf-Rayct companion (van Kerkwijk et al. 1992). As demonstrated in LP96, 
the characteristics of the X-ray line spectrum, including the shape and magnitude of the emission measure 
distribution conform to what one would expect according to a simple model of a photoionized stellar wind 
in an X-ray binary system. 

The purpose of this paper is to extend and refine the analyses of LP96 and N94, and to demonstrate 
that the X-ray spectrum of Vela X-1 can be interpreted in terms of an explicit stellar wind model. Our 
approach improves upon the method of LP96 in several ways. First we take explicit account of wind 
acceleration according to a modified Castor, Abbott, & Klein (1975, hereafter CAK) wind model. We 
model the orbital phase dependence of the apparent differential emission measure, taking advantage of 
the fact that the system parameters (binary separation, companion radius, inclination angle, ephemeris, 
etc.) are much better constrained in Vela X-1 than in Cyg X-3. The LP96 method of fitting the data with 
model ionic spectra, each of which is assumed to be at a single temperature, is adopted as an initial step 
in the course of our analysis. But we also analyze the data with a more elaborate model that accounts 
for temperature variations in each ionization zone. Using the Hatchett & McCray (1977, hereafter HM77) 
picture of an X-ray source immersed in a stellar wind, coupled with a generalized CAK wind, we show 
that the shape and magnitude of the differential emission measure (DEM) distribution in HMXBs depend 
sensitively on stellar wind parameters. Prom a grid of model DEM distributions, theoretical spectra are 
generated and compared directly to the spectral data, leading to the determination of the rate at which 
highly ionized material is accelerated and the corresponding mass loss rate. 

This paper is organized as follows. In §2, we describe the data reduction performed, the details of our 
explicit spectral model, and the derivation of the empirical DEM distribution for the ionized stellar wind. 
In §3, we discuss physical models for the stellar wind, and our derivation of wind parameters. In §4, we 
present a quantitative analysis of the K-shell fluorescent lines. Finally, in §5, we discuss the implications 
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of our results in terms of previous models for the Vela X-1 system, and show that all of the data can be 
reconciled through a model of an inhomogeneous wind were most of the mass is contained in relatively dense 
cool clumps, while most of the volume is occupied by a lower density, more highly ionized tenuous medium. 

2. Spectral Analysis 

2.1. Data Reduction 

Vela X-1 was observed with ASCA on three occasions: twice in 1993 and once in 1995. The first 
set of observations in 1993 was performed during the performance verification phase, and covered the 
entire transition from ingress to egress (—0.13 < (p < +0.16). The observation in 1995 covered only part 
of the eclipse phase (—0.05 < (j) < +0.10). For the purpose of our analysis, we only use data obtained 
in 1993. We separate the observations into three phase ranges: pre-eclipse (—0.13 < (f) < —0.10), eclipse 
(—0.10 < (p < +0.10); and post-eclipse (+0.12 < (f> < +0.16). For the eclipse phase, we use data from SISO 
and SISl. For the pre- and post-eclipse phases, which we analyze principally to derive the intrinsic X-ray 
continuum spectrum of the compact source, we only use data from GIS2. 

The data reduction was performed with FTOOLS v4.0. The data were screened through standard 
criteria: rejection of hot flickering pixels and data contaminated by bright earth, event selections based 
on grade, etc. The observation during eclipse was performed in 4-CCD mode with all chips collecting a 
significant number of counts. In order to maximize statistical quality, source and background counts were 
separately extracted from each individual chip. Data, background, response matrix, and ancillary response 
matrix files were generated for each CCD chip and were combined using standard methods. The resulting 
spectra contained 1.1 x 10^ and 9.0 x 10^ counts from SISO and SISl, respectively, collected during a total 
exposure time of 38 ksec. The net exposme times of the pre- and post-eclipse phases were 6.0 and 8.4 ksec, 
which contained 2.5 x 10^ and 7.4 x 10** counts, respectively. 

2.2. Spectral Model 

Our spectral model for Vela X-1 is comprised of three principal components: (1) continuum radiation 
from the neutron star, which consists of both direct emission from the central source, absorbed through 
the intervening wind, and scattered radiation from free electrons in the extended wind; (2) discrete 
recombination line and radiative recombination continua (RRC) from highly ionized ions; and (3) 
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fluorescent lines from near-neutral ions. Our treatment of each of these components is described in the 
following subsections. 

2.2.1. Continuum Emission 

The spectra of Vela X-1 require two continuum components. The intrinsic continuum radiation, which 
originates near the neutron star, is highly absorbed as it propagates through the surroimding stellar wind. 
A fraction of this radiation is scattered into our line of sight through electron scattering and is less absorbed 
relative to that of the direct continuum. The direct continuum is also required by the eclipse spectrum, 
since near eclipse ingress {(f) ~ —0.10) and egress {(p ~ +0.10), the hue of sight to the compact object is not 
entirely blocked by the companion. The monochromatic photon flux (photon cm~^ s~^ keV~^) of the total 
observed continuum radiation, therefore, has the following functional form. 

F,o„t(£;) = A--*e-'^(^)^-"i;^^ + A'i''-e-(^X"£;,-^^ (1) 

where A^'^^*^ and A*^" are the normalizations (in imits of photon cm^^ s~^ keV^^ at 1 keV) of the scattered 
and direct components, respectively. Since electron scattering should not appreciably alter the spectral 
shape of the continuum, the photon indices of the two components are set equal. The absorption cross 
sections are adopted from Morrison & McCammon (1983) and we assume solar abundances for the absorbing 
medium. 

2.2.2. Discrete Recombination Emission 

Our recombination cascade model is an outgrowth of that described in LP96, which was successfully 
applied to the analysis of the ASCA spectrum of Cyg X-3. Several updates have been made, and the 
current version includes line emissitivies of K-shell carbon, nitrogen, oxygen, neon, magnesium, silicon, 
sulfur, argon, calcium, and iron, as well as the iron L species (Fe XVII - XXIV), and the associated RRC. 
We have written an interface to the XSPEC spectral analysis package (Arnaud 1996) that calculates line 
and RRC emissivities for a given ionization stage and electron temperature. Each ion is treated as an 
individual model component in XSPEC, with the temperature and normalization as adjustable parameters. 

Atomic structure and transition rates are calculated with the atomic physics package, HULLAC 
(Hebrew University/Lawrence Livermore Atomic Code) (Klapisch et al. 1977). Ionic spectra are calculated 
for each ion using radiative cascade models that include detailed atomic structure through n = 6 and 
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averaged structure for n = 7 - 10. The n = 10 level is a "super level" that accounts for recombination 
into shells with n > 10. Radiative rates, including El, E2, Ml, and M2 transitions are calculated with 
HULLAC. For the 2El two-photon decay rates of hydrogenic 2si/2 states, we use Parpia & Johnson 
(1982). We assume an dependence in the photoionization cross sections, where threshold cross sections 
and ionization potentials are adopted from Saloman, Hubble, & Scofield (1988). The RRC are calculated 
from the photoionization cross sections using the Milne relation, assuming that the recombining electrons 
are Maxwcllian distributed. The line and RRC emissivities are self-consistently calculated for each ionic 
species to produce the correct line- RRC ratios. 

The specific line power 5„; (T) for a radiative transition connecting an upper level u to a lower level I is 
given by 5'„;(T) = rjui oib,r{T), where r}ui is the fraction of recombinations onto the relevant charge state 
that produce the transition u ^ I, and ana is the total radiative recombination (RR) rate coefficient for 
that charge state. The dimensions of Sui, like aim, are cm^ s~^. We calculate Sui by solving the set of 
coupled rate equations, assuming that only the ground state of the recombining ion is populated. This is an 
excellent approximation for K-shell ions at most astrophysical densities. For the recombined ion, however, 
the level populations (dimensionless) are computed explicitly. Thus Sui is simply the photon line power 
per ground state of the recombining ion. For a transition with energy Eui , Sui is related to the conventional 
recombination line power Pui (erg cm'^ s~^) by Pui{T) ~ {uh /ne)Az ,fi+i{T)Sui{T)Eui, where Az is the 
elemental abundance and /j+i is the ionic fraction of the recombining ion. The quantities nn and Ug are 
the hydrogen and electron densities, respectively. 

Operationally, Sui is found directly from the inversion of the rate matrix according to 
Sui{T) = Nu{T)Aui/ne, where Aui is the radiative decay rate corresponding to the transition 
u ^ I. The rate matrix is solved over a grid of temperatures, and, for each X-ray transition in the model, 
Sui is fit to a power law; <S'„;(T) = C„; T"'''"', which provides an adequate description for our purposes. 
Values of jui fie in the approximate range 0.6 - 0.8. 

In constructing the spectral model, a number of physically plausible approximations are invoked. 
First of all, we do not include line flux contributions from coUisional excitations since we expect the local 
temperatures to be much lower than the ionization potentials {kT/x ^1). CoUisional transfers in the 
He-like ions are also not included. This approximation is valid in the low density limit (n ^ 10^° cm~^) 
and, in any case, the ASCA CCDs do not have the resolving power to separate these lines to observe 
variations of line ratios in He- like ions. We also do not take into account inner-shell ionization of Li-like 
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ions, which can lead to the production of the forbidden Une in He-Uke ions. 

Although they arc available in our spectral model, we do not include the forest of iron L transitions 
in our analysis of the Vela X-1 spectrum. For X-ray photoionized plasmas, the line power of individual 
iron L lines is low compared to those of H-like and He-like ions of lower Z elements which exist at similar 
ionization parameters (Kallman et al. 1996). In addition, at ASCA SIS resolution, the iron L lines are not 
resolvable - they produce a more-or-less smooth continuum component between 0.7 and 2.0 keV, which 
provides only a minor perturbation to the true continuum in the same energy range. As we show later, the 
omission of these features has a negligible effect on our derived wind for the Vela X-1 system. 

As indicated earlier, the line power of a particular transition produced by recombination cascade 
depends on temperature through the temperature dependence of the recombination coefficient. In principle, 
the temperature appropriate to the relevant charge state can be inferred from the shape of the associated 
RRC feature, as was done using ASCA data for Cyg X-3 (LP96) and 4U 1626-67 (Angehni et al. 1995). 
However, for the Vela X-1 data set, the RRC are blended with other features, and the spectra are 
of somewhat lower statistical quality, so this is not practical in this case. Instead, we must assume a 
temperature for each ionic species. A natural choice for the temperature corresponds to that where line 
powers of the brightest lines for that ion have their maximum. A large fraction of the total ionic line 
emission originates near this temperature. Therefore, for each ion, we fix the temperature at the value 
where the product fi+i{T)Sui{T) peaks. In photoionized plasmas, however, the charge state distribution of 
the gas is determined not by the temperature, but by the ionization parameter, ^ = Lx/{npr'^), where 
is the X-ray luminosity, Hp is the proton number density, and Vx is the distance from the ionizing source 
(Tarter, Tucker, & Salpeter 1969). We use XSTAR to calculate fractional ionic abundances, fi+i, and 
temperatures, T, as functions of ^. For each ion in the model, we then use the T(^) relationship, and fix 
the electron temperature according to the value of ^ for which fi+i{C)Sui{C) attains its maximum. For each 
charge state i, we define the ionization parameter of formation (^form) to be that value of ^ corresponding 
to the maximum of fi+iSui- We assign the temperature of formation according to Tform = T'(Cform)- 

XSTAR requires an input shape for the ionizing continuum X-ray spectrum. At energies below 10 keV, 
this is constrained by the shape of the continuum visible in the ASCA spectrum itself, but the ionization 
and temperature structure is also sensitive to the shape of the spectrum at higher energies, outside of the 
ASCA bandpass. We assume the spectrum determined for Vela X-1 using HEAO 1 data by White, Swank, 
& Holt (1983), which involves a broken power law with an exponential cutofi^ at 20 keV with an e- folding 
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energy of 16 kcV. In Table 1, we list the ^torm and their respective Tform values for the H-like and He- like 
species of oxygen, neon, magnesium, silicon, sulphur, and iron, given this assumed ionizing spectrum. 

EDITOR: PLACE TABLE |l| HERE. 

2.2.3. Discrete Fluorescent Emission 

We do not explicitly model the atomic physics of fluorescent line production, but rather, represent 
them as isolated gaussian components. Apart from the iron K line complex in which the characteristic line 
width can be determined from the data, the mid-Z fluorescent line widths are fixed at cr = 50 eV. Their 
centroid energies and normalizations are left as free parameters. The lines widths were chosen to represent 
typical spreads in fluorescent line energies for the low charge states of each element. In fact, the line fluxes 
are not sensitive to our choice of a as long as they lie within realistic values. 

2.3. Results of Spectral Fitting 

Spectral fits were performed to the SISO and SISl data during eclipse phase with XSPEC. We 
have modified the code to incorporate the specific spectral model described in §2.2.2. Parameters for 
the fit include the normalizations, absorbing column densities, and power law indices for both the direct 
and scattered continuum components, the emission measures (proportional to the normalization of the 
recombination line spectrum and associated RRC) of the H-like and He-like charge states of oxygen, neon, 
magnesium, silicon, sulfur, argon, calcium, and iron and the normalizations and line centroid energies of the 
K-shell fluorescent features associated with magnesium, silicon, sulfur, argon, calcium, iron, and possibly 
nickel. Most of these are left "free" to be constrained by the fit, however, we forced the power law index of 
the scattered component to agree with that of the direct component, and applied the attenuation associated 
with the derived absorbing column density of the scattered component of the line and RRC radiation as 
well. The latter is justified if both the scattered continuum and all of the discrete features originate in the 
same extended stellar wind, as we have assumed. 

Despite the large number of free parameters, there are still many (262) independent degrees of freedom. 
Nevertheless, we find an excellent fit with this model with Xr = 0.88. The data and best-fit model are 
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shown in Figure 1, along with the residuals. Note that all of the subtle discrete structure in the spectrum 
is accurately reproduced. 

EDITOR: PLACE FIGURE |l| HERE. 

EDITOR: PLACE TABLE | HERE. 

The derived line fluxes and centroid energies for the fluorescent features are listed in Table 2. The 
detection of many of these features was previously reported by Nagase et al. (1997). The centroid energies 
suggest that most of this emission arises in material that is partially ionized, as indicated by the likely 
ion identifications listed in the second column of the table. Partial ionization is indeed expected for even 
dense components of the wind, given the intense ultraviolet field of the companion star which irradiates this 
material. The reality of each fiuorcsccnt feature can be assessed by the change in ^ which results when it is 
omitted from the fit, as listed in the fourth column of Table 2. Except for the nickel line at 7.52 keV which 
is required at a 91.7% confidence level, all of the features are highly significant ( > 99.9% confidence level). 

Due to significant line blending at ASCA resolution, there is ambiguity in the derived fluxes of the 
fluorescent lines. While the sulfur, argon, and calcium fluorescent features are well-separated from the 
bright recombination lines and RRC from the more highly ionized material, the magnesium and silicon lines 
are blended. For example, the magnesium line at 1.30 keV lies close to Mg XI Ka (1.34 keV) and the Ne X 
RRC edge (1.36 keV). If we allow the temperatures of Ne X and Mg XI to be additional free parameters 
in the spectral fit, the necessity for the 1.30 keV feature can be removed by lowering the temperature 
of Ne X, which has the effect of increasing the contrast in the RRC with respect to the corresponding 
lines. However, the required temperatures we obtain from this procedure (1.1 eV for Ne X and 5.3 eV for 
Mg XI) are factors of ^ 100 and ~ 4 lower than what would be expected from photoelectric heating in 
photoionization equilibrium at this level of ionization. Since these temperatures are unphysical, we choose 
to fix the temperatures at their respective temperatures of formation to derive the fluorescent line fluxes 
and centroid energies. 

Similar spectral fitting procedures were performed for the GIS2 data during the pre-eclipse and 
post-eclipse phases. Our goal, in these cases, was merely to characterize the continuum so as to better 
constrain the ionizing spectrum and luminosity. Thus, although we include contributions from discrete 
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emission in the spectral fit, we do not attempt to quantify the detailed characteristics of the line spectrum, 
which are poorly determined anyway, given the lower spectral resolution of the GIS. The best- fit continuum 
parameters for the three phases are listed in Table 3. The derived power law indices agree well with those 
determined by previous experiments. Note that the absorbing column density inferred for the scattered 
component of (5 — 10) x 10^^ cm~^ is roughly consistent with the interstellar reddening to this source 
observed in the optical and ultraviolet (Nandy, Napier, & Thompson 1975; Conti 1978). Adopting a 
distance to Vela X-1 of 1.9 kpc (Sadakane et al. 1985), we derive an average X-ray luminosity, which we 
define as the contimmm luminosity above E = 1 Ryd, of = 4.5 x 10'^^ erg s^-^. 



EDITOR: PLACE TABLE | HERE. 



2.4. Derivation of the Empirical DEM 

Since radiative recombination is a two-body process, the emissivity of each recombination line and 
RRC feature is proportional to density squared: juiiO — "-e Pul{£,)^ where P^i is the line power, defined 
earlier in §2.2.2. If we define the emission measure to be EM — J dV, the line luminosity produced 
at the source over a narrow range of ^ (or logf) — the differential line luminosity — can be expressed 
as dLui ~ PuiiO [d{EAI) /d\og£] dlog^. In HMXBs, and probably most accretion powered sources, the 
observed spectrum is likely to be a superposition of individual spectra, each of which can be characterized 
by a single ionization parameter, provided that radiative transfer does not deform or attenuate the ionizing 
spectrum. Therefore, the line luminosity that one would infer from a measurement of the line flux is given 

by 



d{EM) 



, PuiiO- (2) 
dlog^ J 

The bracketed quantity is known as the differential emission measure (DEM) distribution. For each value 
of ^, the DEM acts as a weighting factor that determines the contributions of each ionization zone to the 
total line flux. Spectroscopic analysis in this context is analogous to the use of a temperature-dependent 
DEM distribution in interpreting spectra from plasmas in coronal equilibrium. DEM analyses are used 
widely for modeling the structure of stellar coronae, for example, and have proven to be an effective means 
by which to extract their global properties (Jordan 1975). 
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If we insert the expression for P„; in terms of Sui (§2.2.2) into equation (2), we have 



Lui = Eui I dlogC 



SuM), (3) 



dlog^ 

where d{EMi^i) ~ n^rii^i dV ^ and where we have used n^+i [nn / ne) fi+iAzn^. Operationally, the 
set EMi+i is determined from fitting to the Vela X-1 spectrum. No constraints are applied to the 
chemical abundances or to the charge state distribution. These factors, as well as the electron density, 
are subsumed by EMi^i. It is important to note, however, that each recombination line has a definite 
temperature-dependent relationship to every other line, as well as the RRC, emitted by that ion, and that 
this is fully accounted for in the model. 

In determining the emission measures from the spectral fit, we make a simple approximation that the 
line power of the brightest line in each ion is a flat function of logi^, over a given range Alog^. For most 
H-like and He- like ions, Alog^ ^ 1. We plot the derived empirical DEM distribution as determined from 
the H-like and He-like lines and RRC of oxygen, neon, magnesium, silicon, sulfur, and iron. 

EDITOR: PLACE FIGURE | HERE. 

For a given line flux, the derived EMs are relatively weak functions of the assumed temperature; 
EMi^i (X T'^ , where 7 ~ 0.7. Nevertheless, there is some arbitrariness in the construction of the empirical 
DEM distribution, and we estimate that our derived DEM values may be uncertain by up to a factor of 2 
or 3. 



3. Physical Model of the Ionized Stellar Wind 

The empirical DEM distribution derived in §2.4 is useful for guiding the interpretation of the spectrum 
in terms of realistic physical models of the ionized stellar wind in the system. By constructing this 
distribution, we have effectively removed the "spectroscopic physics" from the problem, and produced a 
curve that can be directly compared to that expected for simple wind models. In this section we explore the 
dependence of the DEM distribution on various wind parameters. This is described in §3.1 - §3.3. Given 
the uncertainties inherent in the derivation of the empirical DEM, however, the final determination of the 
constraints on wind parameters can only come from direct comparison of fully self-consistent wind spectral 
models with the ASCA spectrum itself. This is described in §3.4. 
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3.1. General Characteristics of DEM Distributions for the Ionized Stellar Wind in HMXBs 

We begin by considering the simple case of a hard X-ray point source located near to a companion star 
with a spherically symmetric stellar wind. We note that there are several mechanisms that may distort 
the wind from spherical symmetry: tidal forces of the compact object which may enhance the mass loss 
rate along the binary axis (Friend & Castor 1982; Stevens 1988), suppression of the wind velocity through 
photoionization (Blondin 1994), non-static accretion (Ho & Aarons 1987; Taam & Pryxell 1988), and the 
presence of a bow shock and accretion wake (Fryxell & Taam 1988; Blondin et al. 1990). However, these 
effects, with the exception of the accretion wake, are relatively localized near the compact object. Since 
our purpose here is to calculate the DEM during the eclipse phases (—0.1 ^ (j) ^ +0.1), their effect on the 
DEM is small. Therefore, we neglect all such effects and assume that the density distribution is spherically 
symmetric around the companion star. We also ignore the presence of the extended stellar atmosphere. 
Sato et al. (1986) adopt an atmosphere in which the particle density decays exponentially with distance 
from the stellar surface. This affects the DEM distribution only below log^ ~ 1, which is outside the range 
constrained by our spectral fits. 

We also assume that the radiation field of the compact object is isotropic, and we ignore transient 
effects that may be produced in the wind by X-ray pulsation. The pulse profile of Vela X-1 is complex 
(McClintock et al. 1976), indicating that the emission is not strongly beamed. Recent ST observations 
showed weak (^ 3%) pulsations in Si IV and N V absorption lines, while the X-ray continuum varies by 
~ 50% from the mean intensity (Boroson et al. 1996). Since the gas density in X-ray emission line regions is 
most likely to be lower than in UV emission/absorption regions, we expect the recombination timescales to 
be comparable to the Vela X-1 pulse period, as will be checked later. The wind, therefore, responds to the 
time averaged flux, rather than the instantaneous flux of the compact object. In addition, recombination 
timescales are also much shorter compared to the timescales for individual ions to traverse ionization zones. 
For these reasons, we ignore the dynamical behavior of the ionization structure of the wind, so that it is 
meaningful to speak of local zones of ionization balance fixed in the system. 

In a stellar wind of massive young stars, the material accelerates radially from the stellar surface 
through momentum transfer due to absorption and scattering of UV continuum radiation in the resonance 
lines from low charge states of the abundant elements (Lucy & Solomon, 1970). We assume that the wind 
has an initial velocity, vq, equal to the thermal velocity of the stellar atmosphere (typically 30 km s~^), 
and accelerates until it reaches a terminal velocity, Woo- Generalizing the results of CAK, we can write the 
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velocity profile as, 



R 

v{r) = Wo + I'oo ( 1 - ^ j , (4) 



where i?, is the radius of the companion, r is the distance from the center of the companion, and /? is a 
parameter which determines the shape of the velocity profile. For typical isolated OB stars, it has been 
shown that (3 ^ 0.8 (Friend & Abbott 1986; Pauldrach, Puis, & Kudritzki 1986). 

To derive an expression for the ionization parameter for all points in the wind, we invoke the relation 
for mass conservation in a spherically symmetric flow (M — Anr'^ fimpnp{r)v{r) = constant, where /i is 
the mean atomic weight). The density can then be written in terms of the mass loss rate, the velocity of 
the wind, and the distance from the companion star. Therefore, the ionization parameter near the binary 
system can be recast in terms of position and stellar wind parameters. 



C(r,r,) =4.3x 10 



2 (-^a;)36(Wcx;)8 / »' 



M_7 \r: 



erg cm s ^ , (5) 



where {Lx)3e is the X-ray luminosity (in multiples of 10"^^ erg s~^), (woo)8 is the terminal velocity (in 
multiples of 10* cm s^^), and M-j is the mass loss rate of the companion (in multiples of lO^^M© yr^^). 
We note that ^(r, r^) is insensitive to vq for vq <ti Woo- The key parameters that determine the distribution 
of ^ are the companion radius i?*, binary separation a, the X-ray luminosity L^, the mass loss rate in the 
wind Af, the terminal wind velocity Voc, and the wind velocity profile parameter /3. We adopt the following 
values: companion radius, i?* — 30 Rq, binary separation, a — 53.4 Rq (van Kerkwijk et al. 1995), X-ray 
Imninosity, = 4.5 x 10'^^ erg s~^ (as determined from the pre- and post-eclipse GIS2 spectra), and a 
terminal wind velocity, Voc — 1700 km s^^ (Dupree et al. 1980). These are listed in Table 4 together with 
the derived wind parameters, as described below. 

EDITOR: PLACE TABLE | HERE. 



EDITOR: PLACE FIGURE | HERE. 

The ionization zones of a spherically symmetric wind, ionized by a point source of X-radiation offset 
from the center of the wind, have approximate bispherical symmetry (HM77). The symmetry is exact in 
the case of a constant-velocity wind. An example of the contours of constant ^ for a Vela X-1 wind model is 
shown in Figure 3. For any particular wind model, we calculate ^(r,rx), and construct a theoretical DEM 
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distribution as follows. We adopt a spherical coordinate system, centered on the companion, with a variable 
radial cell dimension increasing with radius like (r — i?*)^, from an initial cell size of Ar = 10^'' Rq. In 
both the polar and azimuthal directions, the cells are divided into 720 bins. The maximum radius of the 
volume is 20 i?,. For each cell, we calculate the electron density from the wind model, and evaluate ^, 
except for those cells lying in the shadow cone, where X-radiation from the compact object is blocked by the 
companion. We separate ^ into logarithmically spaced bins, with each bin having a width of Alog^ — 0.025. 
The DEM contribution is evaluated for each cell, and added to the appropriate ^-bin. Convergence to 
within a few percent at all values of log^ is achieved with this procedure. We refer to the DEM distribution 
calculated in this way as the intrinsic DEM distribution. 

For comparison to the empirical DEM distribution we need to take explicit account of the occultation 
of various parts of the wind by the companion star, which is a function of inclination angle and orbital 
phase. We refer to the phase-dependent, visible portion of the intrinsic DEM as the apparent DEM 
distribution. The phase dependence is shown in Figure 4, where apparent DEMs at three orbital phases 
are plotted. Note that the calculations are carried out for an assumed inclination angle of 74°. This 
permits a partial view of regions of high ^ near the X-ray source, regions that would be entirely occulted 
if the inclination angle were sufficiently close to 90°. Since regions inside the shadow cone and within the 
occulted region must be excluded when calculating the apparent DEM distribution, the total apparent EM 
(integrated apparent DEM distribution) reaches a minimum at cj) — 0.0, when there is minimal (although 
non-vanishing) spatial overlap between these regions. Conversely, the apparent EM reaches its maximum 
value when there is maximal overlap, i.e., when (f) = 0.5. This behavior is illustrated in Figure 4, along with 
the intermediate case, — 0.25. The DEM increases monotonically for all ^ between phases 0.0 and 0.5. 

EDITOR: PLACE FIGURE | HERE. 

3.2. (3 Dependence of the DEM Distribution 
EDITOR: PLACE FIGURE | HERE. 

The /3-dependence of the DEM distribution is illustrated in Figure 5, where we have plotted the 
apparent DEM distribution for the eclipse phase for various values of /?, assuming a mass loss rate, 
terminal velocity, and X-ray luminosity as given in Table 4. For a constant velocity wind ((3 = 0), the 
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DEM peaks sharply at the value of ^ corresponding to the midplane (r = r^) of the binary system — the 
plane perpendicular to the line of centers at the midpoint. We refer to the position of this peak as ^mid- 
The peak owes its existence to a purely geometrical effect; the volume that covers a given range in ^ has 
a global maximum at ^ = ^mid) since the density is a smooth function of position, the DEM peaks at 
that value. The narrow width of the DEM distribution (Alog^ « 0.6) results in an emission line spectrum 
dominated by only a few ions. The f3 = case is the simplest, and does not account for the dynamical 
effects associated with the UV field of the companion. It may, however, have some relevance to the case of 
a "coasting" wind, where the ions responsible for driving the wind have been ionized away by the combined 
UV and X-ray fields (MacGregor & Vitello 1982). 

For an accelerating wind, (/? > 0), the wind velocity near the companion is significantly lower than 
the terminal velocity. For a given mass loss rate, the density in this region is thus higher than that of 
the constant- velocity case. Since, locally, d{EM) oc v~'^, the DEM increases for ^ ^ ^^^^ (the companion 
half-space) as /3 increases. In fact, even for the /? = 0.2 case, which corresponds to extremely rapid 
acceleration, the DEM distribution has started to broaden, with a rapid increase in the DEM near the 
lower edge of the distribution. All traces of a peak in the distribution have vanished for the intermediate 
values of (3. For /3 in the range 0.4-0.7, the distribution resembles a trapezoid, with a width Alog^ near 
1.2. For more slowly accelerating winds (/? ^ 0.8) a peak develops near log^ — 2.0. Comparison with the 
empirical DEM distribution depicted in Figure 2, suggests that the value of /? for the Vela X-1 wind lies in 
the range 0.4 - 0.8. 

Regions with ^ ^ ^mifj (the X-ray source half-space) are nearly unaffected by variations in /3, since, by 
the time the material reaches the X-ray source half-space, it acquires a large fraction of its terminal velocity 
and, therefore, is nearly independent of the value of /3. For higher values of ^, say log^ ^ 4.0, the surfaces 
of constant (, approach perfect spheres concentric with the X-ray source. Moreover, in this idealized smooth 
wind approximation, the density does not vary drastically across a surface of constant ^ in this region(i.e., 
rie is independent of (3). Therefore, the asymptotic form of the intrinsic DEM distribution is given by 
d{EM)/d\og^ = 27rlnl0 [L^Upf/^ ^-3/2. 

Figure 5 also shows that the total EM (the integrated DEM) increases monotonically with (3. For this 
example, EM changes by a factor of 8 in comparing the (3 = 0.0 and /3 = 1.0 distributions. However, for 
the most likely range of values, say f3 = 0.4 - 0.8, EM increases by less than a factor of three. We can use 
this fact to derive a simple way to estimate the mass loss rate, as discussed in the next section. 



3.3. Total Emission Measure 



The total emission measure is simply the volume integral of over the wind. The electron density, rig, 
is related to the mass density, p, by: Ug = Kp/rUp, where k is the mean number of electrons per nucleon. In 
general, k depends on the local value of the ionization parameter. However, for near cosmic abundances, 
free electrons come mainly from hydrogen and helium, so that if these two elements are fully ionized 
throughout the wind, then k is only very weakly dependent on position. In that case, for a spherically 
symmetric smooth wind with a CAK velocity profile, we get 

M^K^ dx 

where x is a dummy variable of integration. 

This is the total emission measure in the wind. To compare with the integral of the empirical DEM 
depicted in Figure 2, we need to convert to an apparent emission measure EM^^pp, which is related to EM 
by a geometric scale factor, / < 1, such that EM^pp = f ■ EM. If we denote the integral by /, equation (6) 
can be inverted for M to give 

1 /2 

M = 2.7 X 10-^ (//)-^/' ( — A ( -S^^l Mq yr-i. (7) 

' Vl700kms-V Vl05^cm-3y °^ ^' 

A lengthy calculation is required to determine / in a general case. For the Vela X-1 system during eclipse, 
/ w 0.30. The integral / depends on /? and the ratio vq/voo- We find that, for typical values of these 
parameters, / falls into the range ^10 - 20. By inspection of Fig. 2, a reasonable estimate for EMg^pp is 
3 X 10^^ cm~^. Therefore, the implied mass loss rate is (1 - 3) x 10"'^ Mq yr~^. 



3.4. Self-Consistent Spectral Fit to the Physical Wind Model 

Having determined approximate values of f3 and M , we can now work backward, and use our physical 
model of the wind to generate an explicit spectral model which can be compared with the ASCA eclipse 
phase data. We assume a specific set of chemical abundances, a detailed model for the charge state 
distribution, and a global model of the wind to determine the density and temperature distribution. For the 
chemical abundances, we take the solar photospheric set from Anders & Grevesse (1989). The ion fractions 
and temperature are calculated for each value of ^ with XSTAR, as described in §2.2.2. The density at 
each point in the model wind follows from specifications of M and f3. The model DEM distributions are 
calculated by averaging the apparent DEMs over the orbital phases —0.1^(1)^0.1, which corresponds to 
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the eclipse phases of Vela X-1, and mimics the time integration of spectra accumulated during an interval 
of changing orbital configuration. 

The 0.0 ~ 5.0 range in log^ is divided into 200 bins of width Alogf = 0.025, as before. At each ^ 
on the grid, T and the set /i+i are found from XSTAR. The line powers of each recombination line in 
the model are evaluated according to our temperature-parameterized line power model. The process is 
repeated for the next value of ^ on the grid, and so forth, until a set of 200 spectral ^-components have 
been calculated. Each DEM model assigns a different weighting to each ^-component; the product of this 
weighting factor and the line power constitutes the integrand of equation (3) . A sum over the log ^-grid 
gives us Lui, as in equation (3). Continuum and fluorescent lines are added to the recombination spectrum. 
The model spectrum is modified by geometrical dilution, and by attenuation by a column density of neutral 
material, then convolved with the SIS instrument response, and compared to the data. Values of are 
determined on a 50 x 70 grid in log A'/ ~ (5 parameter space. The domains of logM [Mq yr^^) and (3 are 
(—6.85, —6.34) and (0.45, 1.14), respectively, with a grid spacing of 0.01 in both dimensions. 

We first perform these comparisons by fixing the continuum fluxes and column densities to their 
best-fit values of the original fit with the recombination cascade model (§2.4). With M and j3 left as free 
parameters, we find that attains its minimum for M ~ 2.5j^Q g x 10^^ Mq yr^^ and /? — 0.651q ;[g, 
where the errors correspond to 90% confidence ranges for two interesting parameters. Although the fit is 
statistically acceptable {xi = 1-27 for 297 degrees of freedom), a better fit is obtained if the column density 
of the scattered component is allowed to vary. The data are then fit with three free parameters: M , /?, and 
7V|j^^'. We find that the best-fit values in this case are M = 2.65l|^Jg x 10"^ Mq yr-\ /3 ^ 0.79i|^;23, and 
^|cat _ 6.8lg'7 X 10^^ cm^^, with xl = 1-01 for 296 degrees of freedom. The errors here correspond to 90% 
confidence ranges for three interesting parameters. The 90% and 99% confidence contours for the two wind 
parameters, M and are plotted in Figure 6, together with some previous determinations using other 
techniques. 



EDITOR: PLACE FIGURE | HERE. 



The best-fit M found here agrees well with the estimate based upon the much simpler approach using 
equation (7). The close agreement is somewhat fortuitous, since our earlier estimate does not take account 
of the details of the DEM distribution. This validates our earlier claim that M can be determined to 
within factors of a few from an estimate of the total emission measure. Our earlier estimate for iJA/aon was 
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3 X 10^^ cm From a detailed spectral analysis, our best-fit value is 2.63 x 10^^ cm 

The best-fit column density is higher than the Galactic value (4 x lO'^^ cm~^) by roughly a factor of 
two. The excess absorption could be ascribed to material local to the Vela X-1 system. Note that /3 and 
^|cat g^j,g ];iig];iiy corrclatcd; since increasing (3 increases the DEM of low-^ material (see Figure 5), a higher 
column density can offset this effect by absorbing the low energy line fluxes. 

As illustrated in Figure 6, our derived range in /3 is consistent with that inferred from UV observations 
(Dupree et al. 1980) and close to the value calculated by Friend & Abbott (1986), while our derived mass 
loss rate is at least a factor of four lower than previously published results, and approximately an order 
of magnitude lower than the midrange of typically quoted values (1-7) x 10~^ Mq yr^^ (Hutchings 
1976; Dupree et al. 1980; Kallman & White 1982; Sadakane et al. 1985; Sato et al. 1986). The best-fit 
model parameters are summarized in Table 4, and the best-fit DEM curve is shown in Figure 7. Also, for 
completeness, the best-fit spectrum and residuals using the three-parameter model are shown in Figure 8. 

EDITOR: PLACE FIGURE HERE. 

EDITOR: PLACE FIGURE | HERE. 

It should be noted, however, that our derived mass loss rate depends on the assumed metal abundances. 
If we assume that the abundances are 1/10 that of solar, the inferred mass loss rate will be higher by a factor 
of \/T0 due to the M^-dependence on the magnitude of the DEM curve. Therefore, M = 2.7 x IQ-'^A/' , 
where is the abundances relative to solar. Similarly, (3 can be decreased if the abundances of low-Z 
elements (O, Ne, and Mg) are increased. The exact dependence of the inferred /3 on Azq is more complex, 
because the /3 dependence on the DEM curve is non- uniform. For example, the DEM at log^ « 2 and 
(3 = 0.8 is a factor of 3 higher than that of /3 = 0.6. Since an increase in DEM can be offset by an 
increase in the abundance of elements for which ^form lies near that value of log^, a large deviation in the 
assumed low-Z metallicity only corresponds to a small change in the derived /3. 

As a cross-check on the wind parameters, we can calculate the implied Thomson scattering depth 
through the wind. Since Thomson scattering is responsible for all of the scattered continuum observed 
during eclipse, the ratio of the normalizations of the scattered and absorbed continuum components 
^scat^^dir gjygg rough cstimatc of the average Thomson optical depth through the surrounding wind. 
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The ratio of the scattered continuum flux to that of the direct continuum is ~ 1.2 x 10~^ during pre-ecHpse 
and ~ 1.6 x 10~^ during post-ecHpse. This is roughly consistent with the angle-averaged electron scattering 
optical depth {{Te)n = 2.3 x 10~^) calculated using stellar wind parameters inferred from the DEM analysis. 
Were M an order of magnitude larger than our best-fit value, the normalization A^'^'^^ would also have been 
an order of magnitude larger. 

In order to justify the elimination of iron L-shell ions from our spectral model, we calculated a model 
iron L-shell spectrum in accordance with the best-fit DEM distribution. Assuming that the iron abundance 
is solar, we find that iron L line and RRC flux is 12% of the total line flux, where the total includes emission 
from all other elements used in the fit. When folded through the SIS response, the iron L-shell spectrum 
resembles a faint continuum, owing to the high line density in L-shell spectra. This introduces a small error 
into the normalizations of the true continuum fiux components, but does not have a significant bearing on 
our results concerning the DEM analysis. The exclusion of iron L-shell ions leads to a great simplification 
in the spectral fit. 

Using the best-fit parameters listed in Table 4, we estimate typical recombination timescales of highly 
ionized ions in the wind. The electron densities in the visible portion of the wind during eclipse are on 
the order of ng 10*~^ cm~'^. For a typical recombination coefficient of aRR ~ 10~^^ cm^ s~^, the 
recombination timescales are longer than 100 s which is comparable to the pulse period of Vela X-1 
(283 s). Since the ionizing emission is not strongly beamed, this justifies our approximation of an isotropic 
radiation field from the compact X-ray source. 

4. K-Shell Fluorescent Lines 

The physical model for the stellar wind that we have used to derive a mass loss rate and velocity profile 
only accounts for highly ionized material (H-like and He- like ions). As indicated, however, our spectral fits 
also require the presence of fiuorescent lines from much lower charge states of many of the same elements. 
The fiuorescent lines have been incorporated into the model on an ad hoc basis, as required by the data. In 

this section, we use the observed fluxes in these lines to infer the conditions in the reprocessing gas which is 
responsible for their emission. 

Fe Ka fluorescence has been previously detected in the Vela X-1 X-ray spectrum with a number of 
earlier observatories (Becker et al. 1978; Ohashi et al. 1984; Lewis et al. 1992), and several candidate sites for 
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its origin have been posited. Ohashi et al. (1984) noted that the iron Une intensity during ecUpse is roughly 
an order of magnitude lower than that at other binary phases, and concluded that a large fraction of this 
emission must be produced within 5 x 10^^ cm of the neutron star. However, in the simple undisturbed 
model of the wind, the local ionization parameter within this region is much too high (log ^ 3) to support 
near-neutral ions. Hence, fluorescent lines must be produced in anomalous, dense structures in this region. 
As shown by Blondin et al. (1990), disruption of the wind due to the presence of the X-ray source produces 
dense filaments in the accretion wake that may be 100 times denser than the imdisturbcd wind. Therefore, 
these filaments have local ionization parameters that are up to two orders of magnitude lower than the 
ambient wind, and arc capable of supporting ions of low charge state. The atmosphere of the companion is 
also a possible site for fluorescent line emission. The atmosphere has densities as high as a few x 10^^ cm~^ 
out to a few tenths of a stellar radius from the photosphere (Friend & Castor 1982). These regions have 
relatively low ionization parameters (Alog^ ^ 1) and enough column density (iVn ^ 10^^ cm~^) to produce 
the observed line fluxes. During eclipse, however, regions near the neutron star, as well as the irradiated 
face of the companion are mostly out of our line of sight. Thus a non- negligible fraction of the observed line 
flux must originate from more extended regions. 

Assuming that the fluorescing medium is optically thin to the ionizing continuum, the luminosity of a 
single K-shell fluorescent line with line energy Ek and fluorescence yield Yk can be approximated by 



where we use the same cutoff power law as in the XSTAR calculations for the shape of the ionizing 
continuum, and a photoionization cross section of the form a{E) = cr{x){x/ , where x is the ionization 
potential. The continuum source has an X-ray luminosity L^- The factor G is a rather complicated function 
that depends on x, -E-c Ef, and F. For a given spectrum, it depends only on x- The final factor in equation 
(8) is the solid-angle averaged photoionization optical depth at threshold. The optical depth here refers to 
paths from the X-ray source to points in the surrounding wind. Explicitly, (t(x))o is given by 



where Nz is the column density of a given element, i.e., a sum over all charge states of a given element. 
The evaluation of equation (9) is simplified by assuming that Nz is angle- independent. The integral then 
becomes Nz^^/^t^, where Ar2/47r is the covering fraction of the fluorescing material with respect to the 
X-ray source. 

We use fluorescent yields from Kaastra & Mewe (1993), and photoionization cross sections from Verner 




(8) 




(9) 
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et al. (1996) to calculate the values of Nz required to produce the observed fluorescent line fluxes. We define 
^equiv ^j^g hydrogen column density of the region which produces the observable fluorescent lines 

(the "equivalent" hydrogen column density), such that Nz = Az N^"" . In evaluating these quantities, we 
assume solar abundances. For the ionizing spectrum of Vela X-1, we find that G — 0.937, which is accurate 
to within 1% for all elements used in our spectral model. 



EDITOR: PLACE FIGURE § HERE. 



We have used equations (8) and (9) and the measured fluorescent line fluxes to infer values of 
N^'^'^{AQ./^t:). The results are plotted in Fi gure 9 versus the atomic number Z of the element involved. 
As can be seen, there is a clear, monotonic dependence on Z. In particular, the derived value of 

is a factor of ^ 20 higher for iron {Z — 26) than it is for magnesium {Z = 12). This obvious 
trend may be a result of a process known as resonant Auger destruction (Ross, Fabian, & Brandt 1996), 
which can be significant for low-Z , L-shell ions in regions of moderate line optical depth. If the optical depth 
is significant, a K-shell fluorescent line emitted by one atom can be resonantly absorbed by another as long 
as there is a vacant hole in the L-shell. Since the autoionization yield for low-Z ions is high in comparison 
to the fluorescence yield, the second excited atom will preferentially Auger decay, thereby destroying the 
original fluorescence photon. Even line optical depths of only a few can almost entirely suppress fluorescent 
Ka emission for such species. The required column densities are typically ~ 10^° cm~^, which is well below 
that inferred from our fluorescent line measurements. 

Resonant Auger destruction is not effective for M-shell ions, where the L-shell is entirely filled and 
resonant K-shell line absorption is prohibited. However, the presence of such very low charge states 
may be suppressed by the intense UV field of the companion star, HD 77581, in the Vela X-1 system 
(Tcff = 26, 000 K; L ^ 5 X 10^ Lq; Hutchings 1974), which can effectively ionize all atoms in the wind with 
ionization potentials less than ~ 50 eV (Kallman & McCray 1982). The low-Z elements: neon, magnesium, 
silicon, and sulfur should exist almost entirely in L-shell species due to this effect, where resonant Auger 
destruction is effective. By contrast, the UV field is not sufficiently intense to fully strip the M-shell of 
higher-Z elements such as calcium and iron. These combined effects at least qualitatively account for 
the lower values of N^^^^ {Ail / in) inferred for magnesium, silicon, and sulfur, and the complete absence 
of fluorescent emission from neon. If this interpretation is correct, then the true column density of cold 
reprocessing gas in the wind is probably closest to that inferred for iron, ^ 10^^ cm^^. 
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As a final comment, note that there are two values plotted for iron in Figure 9, one derived from the Fe 
Ka line, and the other from the Fe K/3 line. The two incommensurate values of N^^^^lACl/iw) arise from 
the apparently anomalous K/3/Ka ratio (0.28lo!o8) compared to the expected value of 0.13 (Kaastra & 
Mewe 1993) under most conditions. A similar anomaly was reported by White, Kallman, & Angelini (1997) 
for the ASCA spectrum of the low-mass X-ray binary 4U 1822 - 37. It should be noted, however, that the 
fluorescent emission line regions are optically thin to both Thomson scattering and Fe K-edge absorption 
as evident from the inferred equivalent hydrogen column density. The spectrum of the primary continuum 
radiation reflected from this medium does not produce an appreciable Fe K-cdge structure and, hence, does 
not contaminate the Fe Kf3 fluorescent line. Therefore, we believe that this anomalous K/3/Ka ratio is real 
and the origin of this effect is still mysterious. 



5. Discussion 

We have shown that the ASCA eclipse phase spectrum of Vela X-1 is quantitatively well described 
by a simple spectral model incorporating direct and scattered continuum radiation from the compact 
X-ray source, hydrogenic and helium-like recombination lines and their associated RRC from a range of 
intermediate-^ elements (oxygen through iron), and fluorescent line emission from near- neutral species of 
many of the same elements (magnesium through iron). The derived ion emission measures for the ionized 
component were assembled into an empirical DEM distribution, which was shown to closely resemble 
predicted DEM distributions from a simple, spherically symmetric CAK model of a stellar wind, ionized 
by a compact X-ray source offset from the center of the wind. We then showed that this same physical 
model of the wind can be used to generate a self-consistent spectral model which also provides a remarkably 
good flt to the ASCA eclipse-phase spectrum. The derived mass- loss rate for the wind from this analysis is 
M = 2.65to 5o ^ 10"^ yr"\ and the derived velocity proflle parameter is /? = 0.79^0:23. 

Although our derived velocity proflle is in good agreement with previous estimates, our derived mass 
loss rate is roughly a factor 10 lower than typical values quoted in the literature for this source. For 
example, Dupree et al. (1980) and Sadakane et al. (1985) found M 1 x 10~^ Mq yr~^ by measuring 
P Cygni line proflles of Si IV and C IV resonance lines during eclipse, and comparing them to theoretical 
predictions. However, it should be noted that the fractional abundance of both Si IV and C IV is quite 
low throughout the wind in their models. Thus their inferred mass loss rates are strong functions of their 
assumed ionization fractions on the "tails" of the charge state distribution curves, which makes them quite 
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sensitive to the assumptions of the model. 

Mass loss estimates have also come from X-ray absorption measurements close to eclipse egress. For 
example, Kallman & White (1982) foimd M = (2 — 4) x 10^*' Mq yr^^, based on cohmm densities inferred 
from Einstein Solid State Spectrometer spectra in the phase interval c6 = 0.19 — 0.22. However, the observed 
phase dependence of the absorbing column density is incompatible with the profile expected from simple 
stellar wind models. In general, the observed column densities are high {Nh ~ 10^'' cm~^) right after eclipse 
egress ((^ ~ 0.1), then decrease rapidly until 0.2, then increase slowly throughout the remainder of the 
orbit. In contrast, smooth stellar wind models predict a steep decrease in absorbing column density only 
within a narrow phase range, 0.10 < (j) < 0.12. Sato et al. (1986) invoked an extended stellar atmosphere 
for the companion star to explain this behavior. In their model, the wind contribution to the absorption is 
small compared to that of the atmosphere for phases 0.10 < cf) < 0.15. However, they still require a mass 
loss rate ~ 10~^ Mq yr~^ in order to get their minimum observed column density, which is 10^^ cm~^, 
significantly above that expected for the interstellar contribution alone. 

Our own best-fit wind model for the ASCA eclipse phase data also requires an absorbing column 
density of ~ 7 x 10^^ cm~^, which is higher than the canonical interstellar value. Since the Vela X-1 system 
is highly inclined {i > 74°; Conti 1978), it is possible that the excess absorption is specific to the equatorial 
plane, where the stellar wind is preferentially enhanced (Blondin et al. 1990). However, it seems unlikely 
that a factor of 10 discrepancy in the mass loss rate can be due to that effect alone, especially during eclipse 
phase, where emission from the vicinity of the equatorial plane is largely suppressed by occupation. 

It is interesting to compare Figure 3 with Figure 9 of Dupree et al. (1980) and Figure 1 of McCray et 
al. (1984), which show values of log^ near the Vela X-1 midplane close to 1.0 rather than 3.0. The higher 
mass loss rates deduced from the UV observations translate directly into larger emission measure at low 
ionization parameters, and are fundamentally incompatible with the observed ASCA eclipse phase spectra 
for smooth stellar wind models. Our calculations show that the DEM distribution at 1.8 < log^ < 2.0, 
for example, would exhibit emission measures that are approximately two orders of magnitude larger than 
derived, which means the expected intensities of O VIII, Mg XI, and Si XIII recombination lines would be 
100 times larger than actually observed. 

This confiict between the previously derived mass loss rates and our own estimate based on the highly 
ionized component of the X-ray spectrum is exacerbated by the presence of fiuorescent emission in the 
ASCA data. Fluorescent emission is the dominant line emission process in lower charge states. However, 
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inspection of Figure 3 indicates that in our smooth wind model there is insufScient low-^ (log^ ^ 1) 
material in the observable portions of the Vela X-1 system to produce these lines. Therefore, irrespective 
of the details of the fluorescent lines, it is clear that their mere existence in the eclipse spectrum signals a 
problem with the simple wind model used to describe the recombination spectrum. 

One way to reconcile the coexistence of highly ionized and near-neutral ions is to consider inhomogeneous 
stellar wind models, where cool, dense clumps of material are embedded in a more tenuous, highly ionized 
wind. Such inhomogeneities are, in fact, expected theoretically for line-driven stellar winds, as a result of 
various instabilities (MacGregor, Hartmann, & Raymond 1979; Lucy & White 1980; Carlberg 1980; Abbott 
1980; Owocki & Rybicki 1984). For example, Rayleigh- Taylor instabilities should produce clumps with sizes 
10^^ cm at intermediate distances from the star (Carlberg 1980). There is also observational evidence for 
the existence of clumps in the wind. Nagase et al. (1986) suggested that partial obscuration, produced by 
clumps with typical column densities ~ 10^^ cm~^, could produce the soft excess observed at particular 
orbital phases in Vela X-1. White, Kallman & Swank (1983) invoked overdense clumps with characteristic 
sizes 10^^ cm to account for observed short-term variability in luminosity and absorbing column density 
for the similar HMXB system 4U 1700 - 37. In that context, it is interesting to note that the second 
ASCA observation of Vela X-1 performed in 1995 exhibited a flare lasting ~ 10^ s, in which the X-ray 
luminosity increased by as a much as a factor of three. It is conceivable that this flare was the result of an 
increase in the accretion rate resulting from a collision of the neutron star with such an overdense region of 
the wind. Since the orbital velocity of the neutron star in Vela X-1 is ~ 300 km s~^, the required size of 
this region is again of order 10^^ cm. 

The presence of cool, low ionization clumps in the wind can be quantitatively constrained by our 
observed fluorescent line intensities. If mass is being lost by the companion star in the form of low ionization 
clumps at a rate Mc, then for a spherically symmetric outflow of these clumps at a constant velocity, v, 
the solid angle subtended by clumps as seen from the compact X-ray source is approximately given by 
fl = '7rMci~l/{mcvRi,), where rc is the typical radius of a clump and rric is a typical mass. For a uniform 
density, the mass of a clump is given by nic = (47r/3)rj?/imj,np, where Up is the proton density within the 
clump, and /x is the mean atomic weight. These last two expressions can be combined to give 

where Nc = «c^c is the column density of a single clump, and we have taken the value for in Table 
4. Equating Nc^/An to the values of A''^''"'^ plotted in Figure 9 for the high-Z fluorescent lines (where 
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resonant Auger destruction is ineffective; see §4), we infer mass loss rates in this clumped component which 
are now quite commensurate with previous estimates. 

The fact that the fluorescent lines are produced in near- neutral species imposes an upper limit, ^max) 
on the ionization parameter in the clumps, which, for a given distance, implies a lower limit on the clump 
particle density: np > 3 x 10^^^~^(r2:/a)~^ cm~^, where rx is the distance from the X-ray source and a is 
the binary separation. For a clump column density of ~ 10^^ cm~^, the implied value of fc is ^ 10^^ cm, 
in good agreement with theoretical expectations. In addition, at the required ^max, the gas temperature 
within the clump is few eV, which means that for the inferred values of rip, the clumps are in near 
pressure equilibrium with the lower-density, higher-temperature wind. 

As demonstrated, the clumped component of the wind accounts naturally for the observed fluorescent 

lines and for the variability in the luminosity and the inferred absorbing column density. Since it enhances 
the column density in low ionization material, it also helps to explain the strength of the low charge state 
absorption lines (C IV and Si IV) measured in the UV. A concern might be raised about the apparent 
smoothness of the observed P Cygni profiles in these lines: since individual clumps should have well-deflned 
velocities, the lines might appear "choppy" if only a few clumps contribute to each absorption line at any 
particular time (sometimes called "discrete absorption components"). However, a simple estimate shows 
that the number of clouds per velocity interval is approximately given by: 



which leads to approximately one clump per 1-2 km s~^ velocity interval. The lUE high-dispersion 
spectra of Vela X-1 had a velocity resolution 20 km s~^, which suggests that the observed profile should 
indeed be smooth. The fraction of the companion disk which should be covered by the clump ensemble is 
approximately: 



where g is a numerical factor of order 1/10. Thus, the large covering factor with respect to the companion 
appears to be easily attainable in this picture. The presence of the clumped component is also supported by 
the orbital phase variability of the P Cygni lines which show evidence of non-monotonicity in the velocity 
structure of the wind (Kaper, Hammerschlag-Hensberge, & van Loon 1993). 

Given our estimates derived above, the clumped component of the wind dominates the mass loss rate 
from the companion. Nevertheless, the hot ionized component dominates the volume of the wind. Within 
an outer radius comparable to the binary separation, we find that the volume filling factor of the clumps is 




2 



(11) 




(12) 
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only 0.04. Therefore, the presence of the clumped component only mildly perturbs the otherwise smooth 
wind of ionized gas, and does not significantly affect our DEM analysis. 

This interpretation of the Vela X-1 spectrum can be quantitatively tested with much higher precision 
when high resolution X-ray spectra become available from AXAF, XMM, and Astro-E. At the resolution 
expected for the spectrometers on these missions, the iron L spectrum will be especially constraining. The 
iron L ions sample a large range in ^ (log^ = 1.8 — 2.8) which roughly brackets the DEM distribution 
inferred from our ASCA observations. The advantage of using iron L lines is that only a single elemental 
abundance must be assumed in establishing the scale for the DEM distribution. In the case of the ASCA 
data, where only hydrogenic and helium-like lines are available, the actual shape of the inferred DEM 
distribution is affected by assumed abundances. 

Future observations will allow the Ka fluorescent line complexes to be resolved into their individual 
charge state contributions. This may allow us to study the acceleration zone of the wind, the atmosphere 
of the companion, and the formation and evolution of the dense structures that we have postulated here. 

We have shown that the orbital phase variations of the apparent DEM can be dramatic for some 
ranges of ^. Observations made at different orbital phases will provide the information needed to refine 
the simple wind model and to describe the extended components of the wind. Substructures such as the 
accretion wake and the magnetosphere may reveal themselves as perturbations to the DEM distribution 
predicted by extrapolating our model distribution to orbital phases away from eclipse. Line flux and 
temperature measurements extracted from orbital phase-dependent, high-resolution spectra, when coupled 
with velocity information obtained from Doppler shifts and broadening, will provide access to a formidable 
set of diagnostic tools. 

This work has benefited from useful discussions with David Cohen, Chris Mauche. and Andy 
Rasmussen. We have made use of data obtained through the High Energy Astrophysics Science Archive 
Research Center Online Service, provided by the NASA/Goddard Space Flight Center. We thank the 
ASCA team and the HEASARC Online Service for their support. This work was supported under a 
NASA Long Term Space Astrophysics Program grant (NAG5-3541). D. A. L. was supported in part by a 
NASA Long Term Space Astrophysics Program grant (S-92654-F). Work at LLNL was performed under 
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Fig. 1. The ^S'C^ SISO data and the best-fit recombination cascade model. The strongest recombination 
and fluorescent lines are labeled. SISl data was also used in the analysis but are not shown for clarity. 

Fig. 2. The derived emission measure distribution based on the best fit to the ASCA SIS data. Emission 
measures from Ar and Ca are not shown in the figure. The error bars correspond to 90% confidence ranges 
for one interesting parameter. 

Fig. 3. — Ionization contours for the parameters shown in Table 4. The contours represent surfaces of 
constant ionization parameter. Axisymmetric surfaces of revolution can be derived by rotating the figure 
about the line of centers. The numbers labeled are log^. The x at (0,0) marks the the position of the 
neutron star, and the shadow cone corresponds to region in which the compact X-ray source is occulted by 
the companion. The coordinate axes are scaled in units of solar radii. 

Fig. 4. — The apparent DEM at three different orbital phases. The top curve is for phase, (f> = 0.5, which 
corresponds to the intrinsic DEM of the system. The middle and bottom curves are the apparent DEM 
curves at quadrature and eclipse of the X-ray source, respectively. The apparent emission measure during 
eclipse is 30% of the total emission measure. 

Fig. 5. — The apparent DEM cmves for = 0.0 to 1.0 in steps of 0.2, averaged over the Vela X-1 eclipse 
phase. All other parameters are shown in Table 4. 

Fig. 6. — The 90% and 99% contour levels for a smooth wind model. Also shown are two sets of wind 
parameters derived using different techniques. 

Fig. 7. — The model DEM that gives the best-fit to the SIS data. The parameters are shown in Table 4. 

Fig. 8. — A comparison of the data with the model generated with the DEM curve shown in Figure 7. The 
resulting Xr is 1-01 for 296 degrees of freedom. SISl data is not shown for clarity. 

Fig. 9. The product of the equivalent hydrogen column density and the covering fraction versus the 
atomic number of the detected fluorescent lines derived from the measured intensities shown in Table 2. The 
two points at Z = 26 correspond to the colmnn densities derived from iron Ka (lower point) and K/3 (higher 
point). A monotonic dependence on Z is observed. 
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Table 1. Energies and Ionization Parameters of Formation for Bright Recombination Lines 



Ion 


Line energy (keV) 


Cform 


Tform (eV) 


VII 


0.562 


1.49 


7 


VIII 


0.654 


1.92 


18 


Ne IX 


0.910 


1.76 


12 


Ne X 


1.022 


2.30 


97 


Mg XI 


1.337 


1.98 


21 


Mg XII 


1.473 


2.58 


123 


Si XIII 


1.847 


2.15 


75 


Si XIV 


2.005 


2.73 


145 


S XV 


2.441 


2.39 


105 


S XVI 


2.622 


2.90 


183 


Fe XXV 


6.667 


3.11 


250 


Fe XXVI 


6.966 


3.66 


650 



"^Centroid energy of the line complex 
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Table 2. Measured Fluorescent Line Energies and Intensities 



Line Energy (keV) Average Ionization Stage Flux (photon cm ^ s ^) Ax^ 



1 "^n+oo^ 

1.0U_Q Q3 


Mg II- VIII 


i2.2toi) X 10-5 


15 


-1 77+0.02 
' ' -0.03 


Sil-V 


(2.01°:^ X 10-5 


21 


f) Qo+0.03 
^•"J^-0.02 


S I-IV 


(2.6l°:^) X 10-5 


33 


^•-'■^-0.02 


Ar II-V 


(2.3t° j) X 10-5 


26 


o 7/1+0.05 
'^-0.06 


Ca I- VIII 


(1.5t0:5) X 10-5 


11 


6 4^^+0 01 


Fc IV-V (Ka) 


(2.51°:^) X 10-' 


242 


7 -1 O+0.04 
' --^^-0.05 


Fc II-V (K/3) 


i7.2t\-l) X 10-5 


31 


7 ro+o.oe 


Ni I- VIII 


{2.2tl:l) X 10-5 


3 



Note. — All errors correspond to 90% confidence ranges for one interesting parameter. 
"^The range in ionization stages are based on the uncertainties in the line energies. 
''The widths of the lines are fixed at a = 50 eV. 

'^The increase in Ax^ when the line is excluded from the spectral model. 
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Tabk^ 3. Mc^asurc'ci Continuum Paranic'tc^rs 



Scattered Direct 



Orbital Phase Nh'^ A" Nh^ T A" 



Pre-eclipse 1.4 1.7 2.4 x 10"^ 58.2 1.7 7.2 x 10"^ 

Eclipse 0.5 1.7 9.0 x 10"* 23.0 1.7 9.1 x lO'^ 

Post-eclipse 1.9 1.7 3.1 x lO'^ 18.5 1.7 5.8 x 10"! 



"^in units of 10^^ cm~^ 

''Fixed to the corresponding absorbed powerlaw index 
'^photon cm"^ s~^ kev~^ at 1 keV 
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Table 4. Adopted and Derived Parameters for the Vela X-1 System 



Parameter 


Value 


Reference 




26,000 K 


Conti 1978 


a 


53.4 Rq 


van Kerkwijk et al. 1995 


Ft* 


30.0 Rq 


van Kerkwijk et al. 1995 


i 


74° 


Conti 1978 


Ec 


20 keV 


White, Swank, & Holt 1983 


Ej 


16 keV 


White, Swank, & Holt 1983 


Lx 


4.5 X 10^^ erg s"^ 


This paper 




0.79 


This paper 




1700 km s-i 


Dupree et al. 1980 


M 


2.7 X lO-'^M© yr-i 


This paper 



-37- 



0.08 




I I I I I I I I I I I I I 

0.5 1 2 5 10 

Photon Energy (keV) 



30 - 



CO 

I 



a 



b 20 - 



o 



10 - 



- 



n 1 1 r 



o 



00 



I I 




I— I tr" 



> 
CO 



CO 



J L 



J I I L. 



J I I L 



Log f 



- 39 - 




-100 -50 



50 

R(R0) 



100 150 



-42- 



-5 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — |- 

Smooth Wind Model 



-5.5 - 





S -6 h 



Sato et al. (1986) 
(X— ray Attenuation) 



* Dupree et al. (1980) 
(UV P Cygni Profile) 



lao 

O 

I — I 

-6.5 I- 




(X-ray DEM) 



-7 



J I I I I I I I L 







J L 



0.5 



1.5 



-43- 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

30 - = 4.5 X 10^^ erg s"^ 

'b ■ M = 2.65 X 10~^ Mq yr"^ 
o 

12 " v_ = 1700 km 



20 - p = 0.79 




J I I I I I I I I I I I I I I L 



12 3 4 

Log f 



-45- 



10 r 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



C\2 

I 



o 

CV2 
CV2 
O 



< 



(D 



1 r 



0.1 r 



t 



10 



J I I I I I I I I I I I I I I I I I L 



15 20 

Atomic Number 



25 



30 



